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Исследованы численные методы решения уравнений, возникающих при использовании Г - распределения для интервального 
оценивания характеристик случайных величин. Приведены результаты решения уравнения с помощью метода Ньютона, хорд и, 
наиболее оптимального, дихотомического деления. 


1. Введение 

В настоящее время уделяется большое внимание 
разработке автоматизированных систем управления 
технологическими процессами (АСУТП). Одной из 
главных задач, возникающих в процессе функцио- 
нирования подобных систем, является задача полу- 
чения достаточно точных оценок неизвестных ис- 
тинных характеристик управляемого объекта [1]. Во 
многих практических случаях при определении дан- 
ных характеристик объектов для обеспечения жела- 
емой точности целесообразно определять не точеч- 
ные оценки их неизвестных значений, а некоторый 
интервал, который с заданной вероятностью накры- 
вает данные значения. Обычно при решении задач 
по определению подобных оценок, особый интерес 
представляет не только точность, но и надежность 


оценок при ограниченном числе измерений. В мате- 
матической статистике решение этих задач произво- 
дится путем построения доверительного интервала 
/„, который с заданной вероятностью а накрывает 
оцениваемое значения [1-4]. 

Доверительный интервал І а устанавливает точ- 
ность определения оценок рассматриваемого пара- 
метра, а характеристикой надежности получаемой 
оценки является доверительная вероятность а. Су- 
ществующий способ построения доверительных 
интервалов с помощью заранее составленных и ши- 
роко известных в математической статистике та- 
блиц малопригоден для применения в АСУТП, во- 
первых, из-за ограниченности значений довери- 
тельной вероятности, для которых составлены дан- 
ные таблицы, и во-вторых, их использование связа- 
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но со значительными затратами памяти компьюте- 
ра для их хранения. Эти обстоятельства обуславли- 
вают необходимость использования других спосо- 
бов построения доверительных интервалов. Как по- 
казывает анализ имеющихся в данном случае воз- 
можностей, существуют два способа, позволяющих 
избежать использования таблиц. Первый из них за- 
ключается в том, что устанавливается функцио- 
нальная зависимость величины доверительного ин- 
тервала от нескольких параметров и, прежде всего, 
от объема выборки значений оцениваемой характе- 
ристики, и далее эта зависимость аппроксимирует- 
ся. Другой из этих способов, предлагаемый в дан- 
ной работе, основан на решении уравнения, возни- 
кающего при построении доверительного интерва- 
ла с помощью того или иного известного численно- 
го метода. Решение названного уравнения в про- 
цессе построения доверительного интервала избав- 
ляет от необходимости использовать ранее соста- 
вленные таблицы и позволяет строить доверитель- 
ные интервалы при любом значении доверительной 
вероятности а и любом объеме выборки. 

2. Идейные основы 

и возможности предлагаемого способа 

Как известно [2, 3], для построения доверитель- 
ных интервалов существуют приближённые и точные 
методы. В данном случае рассматривается точный ме- 
тод. Для точного нахождения доверительных интерва- 
лов необходимо знать заранее вид закона распределе- 
ния случайной величины X, тогда как для применения 
приближенных методов это не обязательно. 

Идея точных методов построения доверитель- 
ных интервалов сводится к следующему. Любой до- 
верительный интервал находится из условия, выра- 
жающего вероятность выполнения некоторых не- 
равенств, в которые входит интересующая нас 
оценка а. Закон распределения оценки а в общем 
случае зависит от самих неизвестных параметров 
распределения случайной величины X. Однако 
иногда удается перейти в неравенствах от случай- 
ной оценки а к какой-либо другой функции на- 
блюдаемых значений х ь Х 2 ,...х„ интересующей нас 
величины X, закон распределения которой не зави- 
сит от неизвестных параметров, а зависит только от 
числа наблюдений и и от вида закона распределе- 
ния величины X. Именно таким образом обстоит 
дело при оценивании математических ожиданий 
случайных величин. Рассмотрим данные возмож- 
ности более подробно. 

Задача построения доверительных интервалов 
для оценки математического ожидания т величины 
X при произвольном числе опытов п решена только 
для случая нормально распределенной случайной 
величины X. Например, доказано [5, 6], что при 
нормальном распределении случайной величины X 
случайная величина Т, определяемая равенством 

т _ Ѵярй-ги) ц. 


где т и 2) - оценки математического ожидания т и 
дисперсии а 2 случайной величины X, вычисля- 
емые в соответствии с равенствами 

і " 1 п 

а) т = Ух и б) В = — Ц-УЧх, ~ »'), 
п " п - 1 " 


подчиняется распределению Стьюдента с л-1 сте- 
пенями свободы. Плотность закона распределения 
имеет вид: 


У ,-.(0 = 


УІ(п-Ѵ)л Г 


п - 1 


^2 Л о 


1+ 


п— 1 


( 2 ) 


где - плотность распределения Стьюдента с 
п— 1 степенями свободы; Г(х) - гамма-функция. 

Для иллюстрации сущности предлагаемого ме- 
тода рассмотрим его применение при построении 
доверительных интервалов для параметра т. 

Пусть произведено п независимых опытов над 
случайной величиной X, распределенной по нор- 
мальному закону с неизвестным параметром т. 
Требуется построить доверительный интервал І а , 
соответствующий доверительной вероятности а. 
Естественно взять этот интервал симметричным 
относительно т. Обозначим через е а половину дли- 
ны интервала. Величину е а нужно выбрать так, что- 
бы выполнялось условие 

Р(\т-т\< е а ) = а, (3) 


где Р(.) - вероятность наступления события, ука- 
занного в скобках. 


Перейдем в левой части неравенства от случай- 
ной величины ( т-т ) к случайной величине Т, ра- 
спределенной по закону Стьюдента. С этой целью 
умножим обе части неравенства \т-т\<е а на поло- 


жительную величину 



В результате, учитывая 


равенство (1), получим: 

С 


Р 


Т <8, 


V 



= а. 


Найдем такое число і а =е а 



, при котором 


выполняется равенство 

Р(|Г| <*„) = «. 


Данное равенство является, очевидно, ни чем 
иным, как уравнением относительно неизвестного 
( а . Отсюда, воспользовавшись плотностью распре- 
деления Стьюдента, получаем 

І а 

Д| 7 ’|<с)= } У,-! (О >Ж = а. 

-І а 


Как видно из формулы (2), - четная 

функция. Поэтому, учитывая аддитивность опреде- 
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ленного интеграла можно видеть, что данное ура- 
внение эквивалентно уравнению вида 

І а 

} $._,(/)* -у = 0. 

В результате, учитывая равенство (2), получаем 
следующее уравнение: 



Задав вполне определённые значения довери- 
тельной вероятности а и числа степеней свободы 
л, получим, что неизвестным в уравнении будет яв- 
ляться величина Величина Г а определяет для нор- 
мального закона распределения случайной величи- 
ны X число средних квадратических отклонений, 
которое нужно отложить вправо и влево от центра 
рассеивания для того, чтобы вероятность попада- 
ния оцениваемого значения т в полученном ин- 
тервале была равна а. Определив і т мы найдём по- 
ловину ширины доверительного интервала 1 а и сам 
интервал 



Как видно из (4), определение интересующего 
нас значения !„ сводится к решению трансцендент- 
ного уравнения, вычислить решение которого ана- 
литически невозможно. Поэтому решить данное 
уравнение можно только с помощью того или ино- 
го численного метода. Как известно, в настоящее 
время существует целый ряд численных методов, 
пригодных для использования и в нашем случае. 
Предварительный анализ данных методов и усло- 
вий их использования в нашем случае позволяет 
видеть, что для решения уравнения (4), наиболее 
целесообразно воспользоваться либо методом Нь- 
ютона, либо методом хорд, либо методом дихото- 


мического деления. В табл. 1 приведены формулы, 
реализующие алгоритмы нахождения величины ?„ 
каждым из этих методов. 

В формулах, приведенных в табл. 1, і ак и і ат — 
соответственно значения /„ на ^-ой и (СН)-ой ите- 
рациях. Процесс вычисления продолжается до тех 
пор, пока на каком-то ^-ом этапе ? аД не станет кор- 
нем уравнения (4) или пока разность между 1 аЛ+] и 
І а к не станет меньше заданной погрешности в, т. е. 
пока не будет выполняться неравенство. 

\х к+1 ~Х к \<Е. 

3. Некоторые результаты применения 

методов Ньютона, хорд и дихотомии 

Все рассмотренные выше методы были приме- 
нены для решения ур. (4). Ниже приводятся резуль- 
таты исследований данных методов, полученные 
при и=30 и а=0,95. Для того чтобы выяснить, как 
ведут себя последовательности приближений, по- 
лученные с помощью каждого из рассматриваемых 
методов, условие окончания процесса вычисления 
решений ! аЬ к= 1, 2, 3... задавалось неравенством 

К* + і-^а.* И 0.00001, 

определяющем достаточно высокую точность ре- 
шений. 

Анализ поведения последовательностей при- 
ближений показывает, что все три итерационных 
процесса в асимптотике сходятся к одному и тому 
же значению { а , вычисленному традиционным ме- 
тодом, т. е. с использованием таблиц значений 1 а 
при а=0,95 и «=30. Однако скорости сходимости 
данных последовательностей существенно разли- 
чаются. Так, по методу Ньютона достаточно точное 
приближение к корню получается на пятой итера- 
ции, а по методу хорд и по методу дихотомическо- 
го деления - соответственно на десятой и на пят- 
надцатой итерациях. 

Скорость сходимости последовательности при- 
ближений зависит для методов хорд и дихотомии 


Таблица 1 . Формулы для решения уравнения (4) различными методами 
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от величины интервала, содержащего корень, а для 
метода Ньютона - от начального приближения і„. 
На рис. 1 и 2 приведены реализации методов дихо- 
томии и хорд, отличающиеся друг от друга длиной 
начального интервала. 



Рис. 1. Поведение последовательных приближений по ме- 
тоду дихотомии при различных интервалах 



Рис. 2. Поведение последовательных приближений по ме- 
тоду хорд при различных интервалах 

Как и следовало ожидать, чем меньше длина на- 
чального интервала, включающего в себя корень 
ур. (4), тем быстрее находится приближенное ре- 
шение заданной точности. 



Рис. 3. Поведение последовательных приближений по ме- 
тоду Ньютона 


На рис. 3 изображено поведение последователь- 
ности приближений по методу Ньютона для 4 О =0 и 
4, 0 =1. Видно, что чем ближе начальное приближе- 
ние 4,о к истинному значению 4, тем меньшее чи- 
сло итераций требуется для отыскания корня ура- 
внения (4). 

Заданная нами достаточно высокая точность 
е=0, 00001 позволяет видеть, что, начиная с пятой 
итерации для метода Ньютона и с семнадцатой для 
метода хорд, наблюдаются колебания в окрестности 
корня. Амплитуда разброса для каждого из этих 
двух методов различная: для метода Ньютона - 
4,6- 10 _3 , для метода хорд- 1,8'ІОА Причина разбро- 
са в том, что при вычислениях на ЭВМ неизбежно 
существуют ошибки округления. Теоретически по- 
следовательности приближений по методу Ньюто- 
на и по методу хорд должны сходиться к истинному 
значению корня. Для устранения данных колеба- 
ний необходимо, очевидно, проводить усреднения, 
и делать это тогда, когда приближенные решения 
уравнения оказываются в окрестности корня реша- 
емого уравнения. Без проведения данных усредне- 
ний или принятия, каких либо иных мер, напра- 
вленных на устранение колебаний, методы Ньюто- 
на и хорд не дают точного решения ур. (4). 

Применение усреднений в методах Ньютона и 
хорд позволяют получить решение, но при этом 
увеличиваются сложность методов и увеличивают- 
ся затраты времени, тогда как использование мето- 
да дихотомического деления позволяет получить 
решение ур. (4) при наименьших затратах времени. 
На основе этого можно сделать вывод, что наибо- 
лее эффективным для решения ур. (4) является ме- 
тод дихотомического деления. 

4. Выводы 

Предложенные алгоритмы решения уравнений, 
возникающих при использовании распределения 
Стьюдента, для нахождения доверительных интер- 
валов, позволяют отказаться от использования та- 
блиц. При этом появляется возможность нахожде- 
ния доверительного интервала для любого задан- 
ного значения доверительной вероятности. 

Рассмотренные алгоритмы решения, являются 
относительно простыми и не требуют больших вы- 
числительных затрат, что дает возможность полу- 
чать доверительные интервалы для оценок параме- 
тров в режиме реального времени. Лучшие резуль- 
таты показал алгоритм, основанный на методе ди- 
хотомического деления. Его использование позво- 
ляет производить нахождение доверительных ин- 
тервалов при наименьших затратах времени и ре- 
сурсов компьютера. 

Рассмотренные алгоритмы могут использовать- 
ся в составе специализированного программного 
обеспечения прецизионных и высоконадежных си- 
стем автоматического контроля и управления раз- 
личными технологическими процессами и объек- 
тами. 
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